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In preparation for the MMS mission we ask the question: 
how common are bipolar signatures linked to the presence of 
electron holes along separatrices emanating from reconnec- 
tion regions? To answer this question, we conduct massively 
parallel simulations for realistic conditions and for the hy- 
drogen mass ratio in boxes larger than considered in similar 
previous studies. The magnetic field configuration includes 
both a field reversal and a out of plane guide field, as typ- 
ical of many space situations. The guide field is varied in 
strength from low values (typical of the Earth magnetotail) 
to high values comparable to the in plane reconnecting field 
(as in the magnetopause) . In all cases, along the separatrices 
a strong electron flow is observed, sufficient to lead to the 
onset of streaming instabilities and to form bipolar parallel 
electric field signatures. The presence of bipolar structures 
at all guide fields allows the control of the MMS mission to 
consider the presence of bipolar signatures as a general flag 
of the presence of a nearby reconnection site both in the 
nightside and in the dayside of the magnetosphere. 



1. Introduction 

NASA's Magnetospheric Multi Scale Mission [Sharma 
and Curtis, 2005; Burch and Drake, 2009] will "use four 
identical spacecraft, variably spaced in Earth orbit, to make 
three-dimensional measurements of magnetospheric bound- 
ary regions and examine the process of magnetic reconnec- 
tion" (quotation from the web site, mms.gsfc.nasa.gov). The 
authors are part of one of the theory teams that is entrusted 
with the task of exploring and quantify specifically the con- 
ditions to be expected for the mission. 

One of the most important tasks is to predict, qualita- 
tively and quantitatively, what signatures are expected un- 
der realistic conditions for typical crossings of the satellites 
in regions of ongoing reconnection. We focus here on one 
specific diagnostics that is among those considered candi- 
dates to flag reconnection: the presence of bipolar signa- 
tures in the electric field. No single diagnostics is likely to 
be a viable flag of reconnection when used by itself, but 
in combination with others, the presence of bipolar electric 
structures can be a valuable signal to consider. 

A key aspect of reconnection is the formation of strong 
electron flows. The electrons become decoupled from the 
ions and can reach speeds greatly exceeding the Alfven 
speed, typical of the bulk plasma motion in the vicinity of a 
reconnection site. The relative drift of electrons and ions or 
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between different electron populations can then largely ex- 
ceed the electron thermal speed and streaming instabilities 
can develop and form bipolar structures [Drake et al, 2003; 
Goldman et al, 2008]. 

Electron holes are localised regions of reduced electron 
density that lead to a bipolar electrostatic field. Bipolar 
signatures in the parallel electric field have indeed been ob- 
served in satellite crossings both in the magnetotail [Cat- 
tell et al., 2005; Pickett et al., 2009] and in the magne- 
topause [Retind et al., 2006]. 

The analysis of the signatures observed by Cluster [Cat- 
tell et al, 2005] led to the conclusion that bipolar structures 
tend to be observed predominantly in the edge of the current 
sheets, along the separatrices (i.e. the field lines emerging 
from the center of a reconnection site) . While bipolar signals 
the relative vicinity of reconnection regions, their presence 
is reported also at a distance of up to 50 ion inertial lengths 
along the separatrices [Retind et al, 2006]. 

Simulations have also shown the presence of bipolar struc- 
tures. Drake et al. [2003] considered an initial force free 
equilibrium that included an initially imposed relative drift 
between ions and electrons (of two times the electron ther- 
mal speed) and a guide field of approximately 5 times the 
reconnecting field. The formation of bipolar structures was 
observed in the vicinity of the reconnection region. Initial 
equilibria more closely resembling parameters observed in 
the magnetic tail, including smaller guide fields of 1/2 of 
the reconnecting field, have been reported by [Cattell et al, 
2005] . The presence of bipolar structures was observed along 
the separatrices leading to the reconnection region. How- 
ever, at guide fields smaller than 1/2 of the reconnecting 
fields, the formation of bipolar structures was suppressed. 
Instead, the same authors observe that in the observations 
bipolar structures are present also in cases where the guide 
field is smaller. 

Cattell et al. [2005] suggest that the disagreement was 
due to limitations in the simulations that were conducted 
for a reduced unphysical mass ratio (m;/m e = 100), a very 
large unphysical Alfven speed (c/va = 20) and a small box 
size (25.6 x 12.8d»). The combined use of large electron 
scales and small boxes introduces strong boundary effects: 
the electron motion along the separatrices is blocked by the 
close proximity of the boundary that, when measured in the 
electron scales, is only 128d e away from the x-point. Cattell 
et al. [2005] proposed that future work revisit the issue re- 
ducing the proximity of the boundary. This can be obtained 
by restoring the physical mass ratio (and thereby making 
the simulation box larger when measured in electron scales) 
and by dicrectly increasing the box size. Both operations 
are rather costly, because more resolution is needed to re- 
solve the smaller electron scale and because of the use of a 
larger simulation domain. 

We have recently developed a massively parallel imple- 
mentation of the implicit particle in cell method, the code 
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iPic3D [Markidis et al., 2010]. iPic3D showed excellent scal- 
ing on massively parallel computers (up to but not limited 
to 16,000 processors) and the ability to conduct simulations 
at the physical mass ratio m,i/m e = 1836 with relatively 
large box sizes (20 x IQdi), resolving the electron skin depth, 
d e = c/cjp e [Lapenta et al, 2010]. The simulations require 
a tremendous computational cost, made possible by a tar- 
geted allocation on NASA's supercomputers for the theory 
support of the MMS mission. 

In the present paper, to address the fundamental need to 
allow the unimpeded electron motion along the separatrices 
for long distances, we extend our previous simulations re- 
ported in Lapenta et al. [2010] by further doubling the box 
size (40 x 20di) in each direction while retaining the same 
accuracy. We used also 4 times as many processors, for a 
total of 4096. With the combined use of physical electron 
scales and of larger simulation boxes, the electrons can now 
move freely for nearly 900d e before the boundary becomes 
an issue. 

The results of the study remove the previous inconsis- 
tency of simulations and observations. When the electron 
motion along the separatrices is unimpeded by the bound- 
ary conditions, strong bipolar signatures typical of electron 
holes are present at any guide field, with the same field in- 
tensities as observed in the satellite measurements. 

This is a crucial result for the MMS mission: bipolar 
structures will not be limited to the cases when reconnec- 
tion happens in presence of significant guide fields. Rather 
their presence is a universal signature of reconnection and is 
therefore a precious requisite for identifying a reconnection 
regions. 

2. Conditions of the Guide Field Scaling 
Study 

We consider an initial Harris equilibrium, with B x — 
Bo tanh(j//i) and p — pt + posech 2 (y / L) , where x is along 
the sheared component of the magnetic field (Earth-Sun di- 
rection in the Earth magnetosphere), y in the direction of the 
gradients (north-south in the magnetosphere) and z along 
the current (dawn-dusk in the magnetotail) . Note the dif- 
ference from the more common GSM coordinate system. 

The chosen physical parameters for the present study are 
consistent with the observed properties of the geomagnetic 
tail: electron thermal velocity v t he/c = .045, ratio of elec- 
tron and ion temperatures T;/T e = 5, current sheet thick- 
ness L/di — 0.5 and background pressure Pb/po = .1. The 
mass ratio is physical, mi/m e = 1836 (neglecting oxygen 
and other minority species). These choices correspond also 
to a realistic Alfven speed c/v a = 300 [Birn and Priest, 
2007]. 

All the simulations use a computational box with L x /di — 
40 and L y /di — 20. Periodic boundary conditions are 
used in x and Dirichelet boundary conditions (E t = and 
<5B„ = 0) in y. The z direction is assumed to be ignorable, 
but all fields and the velocities of the particles have also the z 
component making the simulation 2D3V. During the evolu- 
tion the electron and ion scales change dramatically and lo- 
cally. We have recently studied the grid resolution required 
to resolve the local electron skin depth in the regions of in- 
terest where the bipolar structures develop [Lapenta et al., 
2010]. The same identical resolution in space and time is 
used here and the reader is referred to the previous study 
for a detailed discussion. 

We have changed parametrically the guide field: B z /Bo = 
0, .05, .1, .25, .3, 1. Reconnection is initiated by a localized 
perturbation [Lapenta et al, 2010] that allows the system 



to evolve naturally without further boundary driving. The 
simulation is continued while reconnection evolves, and it 
is stopped before the process starts to be affected by the 
boundary conditions and saturates. 

The resulting history of the reconnected flux and of the 
reconnection rate are reported in Fig. 1. As previous studies 
have indicated [Ricci et al, 2004], the reconnection rate re- 
mains strong for all guide fields, with a progressive decrease 
of the reconnection rate with increasing guide fields that be- 
comes significant when the plasma f3 in the guide field drops 
below unity. [Rogers et al, 2001, 2003]. 

Figure 2 and 3 summarize the results of our investiga- 
tion. The parallel electric field (£j| = E • B/B) is shown 
on the left and the corresponding parallel electron speed 
(I7 e || = U e • B/B) is shown on the right for each different 
guide field. Given the different rate of reconnection for the 
different runs, we choose for each run a time correspond- 
ing to the same reconnected flux (indicated in the figure by 
the horizontal line). The choice of the time during the re- 
connection process is not critical. As reported in Lapenta 
et al. [2010], the electron flux and the bipolar signature are 
present with similar properties in all phases, at the peak of 
the reconnection rate (the time chosen in the present study) 
on later when it reaches a slightly lower value. Two movies 
of the typical evolution of one of the runs (with B z /Bo = .2) 
is reported in the supporting material online, for the elec- 
tron velocity and electric field parallel to the local magnetic 
field. 

Two regions of strong electron flow can be identified: a) 
the electron diffusion region located right at the x-point with 
the electron jets ejected from it; b) the incoming separatrices 
where the electrons are sucked into the diffusion region. 

The diffusion region and the outgoing jets are more 
prominent at low guide fields and are characterised by hot 
electrons with local thermal speeds comparable with the flow 
speed. Under these conditions streaming instabilities can- 
not easily develop while other weaker instabilities might be 
present [Goldman et al, 2008]. However, at the x-point and 
its immediate vicinity, the parallel electron flow is predomi- 
nantly in the out of plane direction and its stability cannot 
be studied by the present 2D simulations. 

The present study, instead, focuses on the separatrices 
where the electron flow has a large in-plane component and 
is the origin of the Hall current connected to the separa- 
tion of scales in magnetic reconnection between electron and 
ions. This topic has been intensely investigated in the past 
and is reviewed in Birn and Priest [2007, chapter 3]. As 
the guide field is increased, the pattern of field-aligned elec- 
tron flow transitions from antisymmetric (horizontally and 
vertically) to symmetric [Rogers et al, 2003; Swisdak et al., 
2005]. At low guide fields, U e \\ is antisymmetric. Recall- 
ing that the upper and lower magnetic fields have opposite 
sign, this leads to a flow that is always moving along the 
separatrices inward towards the reconnection region. At 
higher guide fields, this pattern changes and the flow re- 
mains inward-directed only along the upper-left and lower- 
right separatrices (for the present choice of the sigh of the 
guide field). Instead, along the other two, it reverses and be- 
comes outward directed. Additionally, the inward-directed 
flow becomes significantly stronger than the outward flow 
along the other two separatrices. This is the reason why 
smaller boxes can pick up the formation of bipolar signa- 
tures only for strong guide fields: there the flow is stronger 
and the boundary conditions affect the electron physics less. 

In all cases, however, the flow speeds exceed the electron 
thermal speed and the strength increases as the guide field 
is increased. Note the different color scale in Fig. 2 and 
3. The flow speed is so considerably larger in the highest 
guide field cases as to advise the use of a different scale for 
a better representation. 
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Under these conditions, streaming instabilities lead to the 
formation of bipolar structures associated to regions of de- 
pleted electron density, referred to in the literature as elec- 
tron holes [Drake et al, 2003; Goldman et al, 2008]. The left 
column of Fig. 2 and 3 reports the observed bipolar struc- 
tures as a function of the initial guide field. The electron flow 
along the separatrices is stronger and leads to the formation 
of bipolar parallel electric field signature. In all cases, bipo- 
lar structures are present, and their strength is not affected 
significantly by the presence of guide fields. However, the 
predominant strength of the inward flow along the top-left 
and bottom-right separatrices tends to increase the preva- 
lence of bipolar signatures along those flow channels without 
suppressing them altogether along the other two separatri- 
ces. 

3. Discussion 

The results above relate directly to the expected signa- 
tures to be observed by the MMS mission. The primary nov- 
elty resides on the choice of the physical parameters used: 
i) a large computational box, ii) the physical mass ratio for 
hydrogen and iii) a realistic value of the ratio of the Alfven 
speed with the speed of light. These three conditions con- 
stitute a first. Previous studies were limited by computa- 
tional resources to small mass ratios or to unphysically hot 
electrons or to small box sizes. Instead, the results pre- 
sented here remove all these limitations thanks to the use 
of the massively parallel code iPic3D and thanks to a large 
computing time allocation on Pleiades devoted to the MMS 
mission. 

The conclusion of the present study can be summarized 
in two main points. 

First, we observe bipolar structures always, for all guide 
fields, including for zero guide field. This is a key prediction 
of relevance to the MMS mission: bipolar signatures associ- 
ated with electron holes are expected to be a common fea- 
ture of reconnection separatrices, regardless of the intensity 
of the guide magnetic field, for component or antiparallel re- 
connection. This extends previous findings that were able to 
find bipolar structures with guide fields from approximately 
5 times [Drake et al, 2003] to 1/2 times [Cattell et al, 2005] 
of the reconnecting fields. 

Second, the present study does quantitative rather than 
qualitative prediction thanks to the physical mass ratio and 
the realistic choice of other parameters. The strength 
and size of the bipolar structures were studied in detail 
in [Lapenta et al., 2010] for the case of B z /Bo = 0.5 and 
it is found here to be insensitive to the guide field value 
for the range considered.. The predicted values of the elec- 
tron speeds and bipolar signature strengths are of primary 



importance. These can be directly compared with satellite 
observations. The bipolar signatures observed have a phys- 
ical extension of approximately 16ADe,; oc (where the Debye 
length is measured with the local density and temperature of 
the regions where the bipolar signatures are found) and are 
of the order eE\\/(miCU) P i) = 1.5- 10~ 5 . [Lapenta et al, 2010] 
These values are in good agreement with the observations 
reporting electric fields of the order of 10-30mV/m [Cattell 
et al, 2005], corresponding to &E\\/ '(rmcwpi) = (1 — 3)-10 -5 . 

The present study has spanned guide fields typical of 
both the magnetotail and the magnetopause but the lat- 
ter presents also distinct asymmetries in the magnetic field 
and density between the solar and Earth side. This effect 
has been neglected here and future work will need to con- 
sider the effect in preparation of the first phase of the MMS 
mission devoted to the magnetopause. Another important 
limitation is that the present study is 2D and flow insta- 
bilities can only develop along the separatrices. To study 
the close proximity of the x-point, where the flow is in the 
out of plane direction, full 3D runs are the next step. This 
limitation is dictated by the need to use the correct physical 
mass ratio, a realistic Alfven speed, sufficient accuracy and a 
large simulation box. These choices have been so far incom- 
patible with a full 3D simulation. Yet the present runs used 
only 4096 cores of the 91,136 available on Pleiades. Newer 
supercomputers have even more. Our future work will focus 
on full 3D study. 




Figure 1. History of the reconnected flux (normalised 
using the in plane lobe magnetic field Bo and density 
rib, left) and reconnection rate (right) for runs at differ- 
ent guide fields indicated in the legends. The horizontal 
dashed line reports the level of reconnected flux used in 
the analysis below. 
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Figure 2. Left: Parallel electric field, normalized as eE\\/(m,iCuj pi ). Right: Parallel electron velocity normalized as 
U e \\/vthe{x,y) (left). From top to bottom, guide fields: B z /Bo = 0, .05, .1, .25. Superimposed in all plots (in white): 
intersections of the magnetic field surfaces (contours of A z ). The time for each plot corresponds to approximately the same 
level of reconnected flux shown in Fig. 1: the actual time chosen is that of the closest graphics dump from the run. Here 
the closest graphics dump is the same for all runs. 
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Figure 3. Left: parallel electric field, normalized as eE\\/(miCu> P i). Right: Parallel electron velocity normalized as 
U e \\/vthe{x,y) (left). From top to bottom, guide fields: B z /Bq = 1/3, 1. See Fig. 2 for details. 
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Figure Captions 

Figure 1: History of the reconnected flux (normalised 
using the in plane lobe magnetic field B and density 
rib, left) and reconnection rate (right) for runs at differ- 
ent guide fields indicated in the legends. The horizontal 
dashed line reports the level of reconnected flux used in 
the analysis below. 

Figure 2: Left: Parallel electric field, normalized as 
eE\ \/(rrii cw p i ) . Right : Parallel electron velocity normal- 
ized as U e \\/vthe(%, y) (left). From top to bottom, guide 
fields: B z /B — 0, .05, .1, .25. Superimposed in all plots 
(in white): intersections of the magnetic field surfaces 
(contours of A z ). The time for each plot corresponds to 
approximately the same level of reconnected flux shown 
in Fig. 1: the actual time chosen is that of the closest 
graphics dump from the run. Here the closest graphics 
dump is the same for all runs. 

Figure 3: Left: parallel electric field, normalized as 
eE\\/ {miOjjpi). Right: Parallel electron velocity nor- 
malized as U e \\/v t he{x, y) (left). From top to bottom, 
guide fields: B z /B a = 1/3, 1. See Fig. 2 for details. 



